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Аннотация 


Введение. Освоение полярных районов Мирового океана, необходимость решения различных задач, связанных 
с наличием большого числа замерзающих внутренних водоемов, ставят перед наукой новые проблемы. К их 
числу относится проблема изучения поведения ледяного покрова под воздействием на него различного вида 
нагрузок. Большой интерес представляет рассмотрение задач о действии на ледяной покров подвижной нагрузки. 
Подвижная нагрузка моделирует действие на лед движущихся транспортных средств. Однако в работах, посвя- 
щенных вышеуказанным задачам, рассматриваются случаи движения нагрузки по прямолинейной траектории. 
Целью данной работы является разработка метода исследования поведения ледяного покрова под действием 
нагрузки, перемещающейся произвольным образом. 

Материалы и методы. В статье предложен метод решения задачи о действии на ледяной покров водоема ко- 
нечной глубины движущейся по произвольной траектории силы. Задача сводится к решению системы двух диф- 
ференциальных уравнений. Первое из них моделирует поведение ледяного покрова и является уравнением коле- 
баний вязкоупругой пластины. Второе — моделирует поведение жидкости, находящейся в состоянии потенци- 
ального течения, и является уравнением Лапласа. Для решения системы дифференциальных уравнений приме- 
нялись интегральные преобразования по временной и пространственным переменным. Полученное в результате 
решение выражалось через повторный интеграл, для вычисления которого применялись численные методы. 
Результаты исследования. В результате реализации предложенного метода получено решение задачи о движе- 
нии сосредоточенной силы по ледяному покрову по произвольному закону. При этом произведены исследования 
характера поведения перемещений и напряжений в ледяном покрове в зависимости от скорости и ускорения дви- 
жения вертикальной нагрузки, глубины водоема и вязкоупругих свойств льда. Кроме того, рассчитано распреде- 
ление вектора скорости частиц жидкости по глубине водоема. 

Обсуждение и заключение. Предложенный метод является весьма эффективным для решения задач о подвиж- 
ных нагрузках, действующих на ледяной покров водоема конечной глубины. Он позволяет решать задачи о дей- 
ствии нагрузки, движущейся по ледяному покрову по сложной траектории. Полученные результаты могут быть 
использованы для расчета напряжения и перемещений ледового покрова при прокладке ледовых дорог или стро- 
ительстве аэродромов на льду. 
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Введение. Освоение полярных районов Мирового океана и наличие большого количества замерзающих внут- 
ренних водоемов приводят к необходимости изучения полей перемещений и напряжений ледяного покрова, обу- 
словленных действием различного вида нагрузок. Решению этих проблем посвящено большое количество работ 
отечественных и зарубежных ученых. Ранее было установлено, что механические свойства льда зависят от соле- 
ности воды и его температуры. Большое внимание уделялось разработке численных моделей льда, которые до- 
статочно точно отражали взаимодействие льда и идеальной несжимаемой жидкости. В работах [1, 2] для этого 
применялся метод гидродинамики сглаженных частиц, в [3, 4] — метод дискретных элементов. В статье [5] лед 
моделируется упругой пластиной, лежащей на поверхности стратифицированной жидкости. Модели, допускаю- 
щие наличие трещин, рассматривались в работах [6, 7]. Модели льда, усиленного армирующими элементами, 
представлены в работах [8, 9]. 

При этом в некоторых работах ледяной покров рассматривается как упругая пластина, лежащая на поверхно- 
сти водоема [10, 11]. В то же время в [12] на основе приведенных исследований делается вывод о том, что в 
некоторых случаях свойства льда лучше всего описывает реологическая модель Кельвина-Фойгта с одним пара- 
метром (временем затухания). Поэтому многие исследователи при моделировании ледового покрова применяют 
вязкоупругую пластину [13]. В [14] для описания свойств льда применялись нелинейные модели. 
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В некоторых работах рассматривалось действие на ледяной покров подвижной нагрузки. В [15] исследовалось 
действие подвижной нагрузки на ледяной покров в замороженном канале, в [16] рассматривалось действие на 
ледяной покров нагрузки с импульсивным характером движения. Работа [17] посвящена изучению нагрузки, дви- 
жущейся по замерзшему руслу. При этом исследовалось прямолинейное движение нагрузки [18]. Однако в ре- 
альных условиях часто приходится иметь дело с нагрузкой, перемещающейся более сложным образом. Поэтому 
целью данной работы явилась разработка метода решения задач о действии нагрузки, движущейся по ледяному 
покрову по сложной траектории. Это позволит более точно исследовать действие на лед транспортных средств, 
движущихся сложным образом. 

Данная работа является продолжением исследований, связанных с рассмотрением задач о действии подвиж- 
ной нагрузки на различные объекты, результаты которых представлены в работах [19, 20]. 

Материалы и методы. Постановка задачи. Рассматривается водоем конечной глубины с бесконечным ледя- 
ным покровом (бесконечная пластина), который подвержен действию вертикальной силы, движущейся произ- 
вольным образом — импульсивно. Предполагается, что жидкость водоема несжимаема и совершает потенциаль- 
ное движение. 

Задача сводится к системе дифференциальных уравнений [15]: 


а = 22275), 
ДЕ = 0, 
где И’, у, ) — прогиб ледяного покрова; Е и и — соответственно модуль Юнга и коэффициент Пуассона льда; 
"= Е Һ/12(1 – 2) — цилиндрическая жесткость на изгиб; й — толщина ледяного покрова; т, — время релаксации де- 
формаций; Д2 = (022 + 22); А= 22 + 0} + 02; рли рь — соответственно плотность льда и воды; с? = р. 0; К = рь 2/0; 
Б = рур; О(х, у, — действующая на поверхности льда нагрузка; Е(Х, у, 2,  — потенциал скорости. 
При краевых условиях при 2=0 (граница лед-вода): 
О.И = 0,Е. 


На дне водоема при 2 = -Н: 

0,Е = 0. 

Кроме того, предполагалось, что ледяной покров и жидкость в водоеме в начальный момент времени находи- 
лись в состоянии покоя. Нагрузка представляла собой сосредоточенную единичную силу (величиной в один 
Ньютон) О(х, у, И, которая перемещалась произвольно по незамкнутой кривой ү произвольной формы. Полага- 
лось, что перемещение силы задавалось в виде О = 0(5(0)), где 5 — дуговая координата, отсчитываемая от неко- 

> е х= хо (1) 
торой фиксированной точки траектории ү. Траектория движения задавалась параметрически в виде , 
у= уо (1) 
где ѓ — время. 
Подвижная нагрузка аппроксимировалась выражением 


О(х,у,1) = е? 1 в ((5 хо (1) +(у= уо (0). 


где = — числовой параметр. 
После применения интегрального преобразования Фурье по переменным х и у, интегрального преобразова- 
ния Лапласа по 1, были получены формулы для вычисления неизвестных функций Ии Е: 


1 Е о | А 1 
УР = 22-р /4= Ј Р -(а-4)т -(а+а)т 
|17 (820) |2 е о (р (3) =е арат, 


нон 


160 _ е0 |арат, 


то 
1 2 2 
(х) | рет" А а): 
0 


о 
С + а)е “+4 =(й= а) ] арат, 


В? (т) =5? +р?, 8 =жо (т) —х, В = уо (т) у, 
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5 


. 12 тор 
у = [0р9 –4(с2р? +врей(рН))(р* +®) | ,а= 2(с2р рН) 


а= | | 
2(с2р +Бет(рН)) 

Используя известные соотношения из теории тонких пластин и теории потенциального течения идеальной 
жидкости, можно получить соотношения для вычисления перемещений и напряжений в ледовом покрове, а также 
компонент вектора скорости частиц жидкости. 

При вычислении несобственного интеграла с помощью численных методов использовалось приближенное 


о А 
соотношение | 7(р)ар я) 7(р)ар, в котором величина А выбиралась настолько большой, чтобы оценка 
0 0 


со 
ошибки | 7(р)ар не превышала установленного значения. 
А 


Так, для величины прогиба льда 


ро ТА ро 
1 1 1 
И (х,у,) = т [7 (т, р)арат = т | бп | [0 (т, р) арат, 


данная оценка имеет вид: 


Г) 
1 24е? 214,2 
— | | Оо (т, р) арат| < С, 
521 о (т р) е хрү(а) 


12 
(4) = [84 -4(с2242 +влет(АН))(А* +). 
Подобные оценки можно получить и для остальных вычисляемых величин. Эти оценки использовались при 
определении величины А. 
2 
2 АЕ А42 [4=? 


В проведенных расчетах величина А выбиралась такой, чтобы оценка арта) 
ды 


не превышала 


1 ГА о 
—— Оз (т, р)арат| 0,001. 
551,1, Д р) А 


При вычислении повторного интеграла использовались квадратурная формула Симпсона (по переменной т) и 
квадратурная формула Чебышева с равными весами для двух узлов (по переменной р). Аналогичным образом 
вычислялись и остальные величины. 

Результаты исследования. Разработан метод решения задач о действии нагрузки, движущейся по ледяному 
покрову водоема, наполненного идеальной жидкостью, по сложной траектории с переменной скоростью. Приме- 
няя данный метод, проведены расчеты, которые показали степень влияния различных параметров на деформацию 
ледяного покрова. 

Изложенный метод не накладывает ограничений на форму траектории движения сосредоточенной силы. При 
расчетах рассматривался частный случай траектории, состоящей из дуг окружностей (рис. 1). Красной точкой 
указано положение сосредоточенной силы в рассматриваемый момент времени и направление движения силы. 


У, ме 
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0 | 4 | 8 12 | 16° Ям 


Рис. 1. Траектория движения сосредоточенной силы 
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При расчетах принимались следующие значения используемых параметров: толщина ледяного покрова 
р = 0,25 м, модуль Юнга Е = 500 000 000 Н/м?, коэффициент Пуассона льда џ = 1/3, плотность льда р = 900 кг/м?, 
плотность воды р = 1 000 кг/м?, = = 5. Ниже представлены результаты проведенных расчетов. 

На рис. 2 приведено изменение прогиба ледяного покрова при скорости движения силы у = 2,5 м/с, касатель- 
ном ускорении у, = 1 м/с?, глубине водоема Н = 25 м и времени релаксации т, = 1 с. 

Закон движения силы по траектории принимался в виде: 

з= 4101 +а2? +а;і. 
Коэффициенты а, а2, аз подбирались таким образом, чтобы сила, находясь в одной и той же точке траектории, 


имела необходимые скорость и касательное ускорение. 


Рис. 2. Изменение прогиба ледяного покрова 


При других значениях указанных параметров качественный характер распределения прогиба ледяного по- 


крова оставался практически неизменным. 


7, м 


У, м 
Рис. 3. Движение жидкости, обусловленное действием подвижной нагрузки на ледяной покров 


Движение жидкости, обусловленное действием подвижной нагрузки при тех же значениях скорости, каса- 
тельного ускорения движения нагрузки, времени релаксации и глубины водоема, представлено на рис. 3 (изоб- 
ражено распределение вектора скоростей частиц жидкости). 


№, м:107 Й м:107 
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Рис. 4. Изменение величины максимального прогиба ледяного покрова в зависимости: 
а — от величины касательного ускорения; 6 — от глубины водоема 
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Влияние на максимальный прогиб ледяного покрова касательного ускорения движения силы представлено на 
рис. 4 а. В этом случае скорость движения силы была равной у = 17,5 м/с, а время релаксации то = 1 с. 

На рис. 4 6 изображен график зависимости величины максимального прогиба ледяного покрова № от глубины 
водоема Н. Скорость движения нагрузки при этом равнялась у = 17,5 м/с, касательное ускорение у; = 1 м/с? и 
время релаксации т, = 1 с. 


ми ——аАы—ы 


—1,0 } 


14. 


и 
0 10 20 30 Г, м/с 
Рис. 5. Изменение максимального прогиба ледяного покрова в зависимости от скорости движения силы 


Зависимость максимального прогиба ледяного покрова от скорости движения силы представлена графически 
на рис. 5. Глубина водоема при этом принималась равной 25 м, а величина касательного ускорения — уу; = 1 м/с2. 
Сплошной линией изображена зависимость, соответствующая времени релаксации т, = 1 с, пунктирная линия 
соответствует времени релаксации т, = 10 с. 

Обсуждение и заключение. Исследовано влияние глубины водоема на максимальный прогиб льда. Получена 
картина прогиба ледового покрова, обусловленная действием сосредоточенной силы, движущейся по сложной 
траектории с переменной скоростью. Расчеты показали, что с увеличением глубины водоема максимальный про- 
гиб ледяного покрова уменьшается (рис. 2). При этом заметная зависимость прогиба ледяного покрова от глу- 
бины водоема Н имеет место лишь для Н < 25 м. При больших глубинах величина максимальных прогибов ста- 
билизируется около некоторого постоянного значения и практически не изменяется. Таким образом, если 
Н> 25 м, то при расчетах глубину водоемов можно считать бесконечной. 

Увеличение касательного ускорения приводит к увеличению прогиба ледяного покрова. Причем зависимость 
прогиба от касательного ускорения очень близка к линейной зависимости (рис. 4). 

При малых временах релаксации т, скорость движения нагрузки заметно влияет на величину прогиба льда. При 
больших временах влияние скорости движения нагрузки на прогиб ледяного покрова заметно уменьшается (рис. 5). 

Для изучения состояния жидкости водоема определено распределение вектора скорости движения частиц 
жидкости, обусловленное действием подвижной силы по льду (рис. 3). 

Разработанный метод решения задач и результаты, полученные с его помощью, могут быть использованы при 
строительстве ледовых дорог, проектировании и строительстве взлетно-посадочных полос на льду. 
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